Nonlinear effect on quantum control for two-level systems 
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The traditional quantum control theory focuses on linear quantum system. Here we show the 
effect of nonlinearity on quantum control of a two-level system, we find that the nonlinearity can 
change the controllability of quantum system. Furthermore, we demonstrate that the Lyapunov 
control can be used to overcome this uncontrollability induced by the nonlinear effect. 
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Quantum control theory is about the apphcation of 
classical and modern control strategy to quantum sys- 
tems. It has generated increasing interest in the last few 
years due to its potential applications in metrology [U J, 
communications [1, 0] and other technologies [il,!^,^^ 3 5 
as well as its theoretical interest. As the effective com- 
bination of control theory and quantum mechanics, the 
quantum control theory is not trivial for several reasons. 
For classical control, feedback is a key factor in the con- 
trol design, and there has been a strong emphasis on 
robust control of linear systems. Quantum system in 
feedback control, on the other hand, can not usually be 
modeled as linear control systems, except when both the 
system and the controller as well as their interaction are 
linear. In fact for many quantum systems, the nonlin- 
ear effects can not be negligible, and in some cases they 
dominate the dynamics of quantum system. Moreover, 
feedback control requires measurement of an observable 
and returns the measured result as a control back to the 
quantum system. This renders the dynamics of the quan- 
tum system both nonlinear and stochastic p^of. In special 
cases the resulting evolution can be mapped into a lin- 
ear classical system driven by Gaussian noise, and con- 
sequently the optimal control problem can be solved by 
classical control theory. However, most control problems 
for such a quantum system can not be solved in this way. 
Therefore a study on nonlinear effects in quantum control 
theory is highly desired. 

Lyapunov functions have played a significant role in 
control design. Originally used in feedback control to 
analyze the stability of the control system, Lyapunov 
functions have formed the basis for new control de- 
sign. Several papers have be published recently to dis- 
cuss the application of Lyapunov control to quantum 
systems [lil, 0, [H, [H, [ll, 0. Although the basic 
mathematical formulism is well established, many ques- 
tions remain when one uses the Schrodinger equation or 
the master equation to describe the dynamics of quan- 
tum system, for example the nonlinear effect in quantum 
control. 

In this paper, we shall address this issue by analyzing 



a two- level system with nonlinear effect. Before studying 
the nonlinear effect in quantum control on the two-level 
system, we recall that a linear two-level system is control- 
lable by two independent parameters, then we show that 
nonlinear interactions may turn the controllable two-level 
system into uncontrollable one. This nonlinear effect may 
result from feedback control, and the uncontrollability 
can be overcome by Lyapunov control as we shall show. 
Consider a two-level system described by. 



(1) 



where cf^ and cfz are Pauli matrices. This model was 
proposed to describe the tunneling of quantum system in 
a double-well potential. In this model, v is the coupling 
constant of the two wells. R denotes the energy difference 
between the two levels. We first show that this system 
is controllable by manipulating the two independent pa- 
rameters R and v. The controllability requires all initial 
states in the Hilbert space Tig of the system can evolve 
to an arbitrary pure target state. This requirement for 
the initial state can be partially lifted by requiring that 
the Hamiltonian H is unchanged up to R and v under 
the following unitary transformation [18], 
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where < < 27r, < ^ < tt. By unchanged we mean 
H' = H{R' ^v') = F^HF, namely the transformation 
F{6^(j)) changes the control parameter in the Hamilto- 
nian H only. The proof is straightforward. Consider 
the Schrodinger equation ih-^\(p{t)) = H\ip{t))^ where 
\^{t)) is the wavefunction of the system. By the time- 
independent transformation FiO^cj)) , \^{t)) F\(j)'{t)) 
we find in§-^{F\^'{t))) = H{F\^'{t))), and ifi^-^W it)) = 
F^HF\(t)'(t)). Since H{R',v') = F^H{R,v)F, we claim 
that there exists a one-to-one correspondence in sets 
{If (t))} and {\^{t))}. Therefore, if {|f (t))} covers ah 
(pure) states in TYs, {|^(^))} is a convex set of all pos- 
sible (pure) states for the two-level system. This obser- 
vation tells us that if the system initially prepared in 
state |e) can be controlled to evolve to an arbitrary tar- 
get state driven by the Hamiltonian H' = H{R' ^v') = 
F^H{R, v)F, the system is controllable. This exactly 
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FIG. 1: (color online) Accessible states (blue points on the 
Bloch sphere) of a two- level system driven by Hamiltonian 
Eq.([T]) in a time interval t G [0,0.5] for (a) and (b), and 
t G [0,4] for (c) and (d). The parameters R and v both 
range from to 7. The initial state is chosen to be |e). All 
parameters are dimensionless. (a) and (b) [(c) and (d)] are 
the same but show the results from the opposite direction. 
P — {Px,Py,Pz) is the Bloch vector, and ^ = = 0. 
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the case as shown in FigdJ where we plot the accessible 
states represented by the Bloch vector p = (Px^Py^Pz)- 
The Bloch vector is connected to an arbitrary state 
\(p{t)) = a{t)\e) + b{t)\g) of the two- level system through 
p=\^{t)){^{t)\ = ^ + ^y-^withp, = 2|a(t)p-l,p^ = 
a*{t)b{t) + h\t)a(t), and Py = i{a{t)b''{t) - b{t)a*{t)). 
We find from Fig{T] that by varying the parameters R 
and the two-level system indeed can evolve to an ar- 
bitrary target pure state, provided the evolution time is 
long enough and there is a wide range of parameters R 
and V to manipulate. It is worth addressing that the 
Hamiltonian in Eq.([T]) becomes H' = ^(7z^\(^x^\(^y 
after the unitary transformation F{0^(j))^ where R' (or 
v') is a function of R^v^O^ and (j). u' is in general not 
zero. At first sight, H does not satisfy the condition 
H{R' ^v') = F^HF, then the two-level system driven by 
Hamiltonian Eq.([T]) is uncontrollable. This is not the 
case, however, because there are only two independent 
parameters {R and v) in the Hamiltonian, hence u' is not 
independent and may be treated as a constant. There- 
fore the term with in plays no role in the control 
on the two- level system[l9|, [20|. 

Now we study the effect of nonlinearity on the con- 
trollability of the system. For this goal, we consider a 
nonlinear two-level model. 
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where \'tp) = \tp{t)) 



C 
2< 



(3) 



( a{t) 



and the parameter C 



characterizes the nonlinear interaction strength, and the 
other parameters have the same notations as in Eq.([T]). 
This model can be used to describe the tunneling of Bose- 
Einstein condensates in a double-well potential and was 
widely used to study the self-trapping and tunneling in 
those systems. 



FIG. 2: (color online) Effect of nonlinearity on the accessible 
states (blue points on the Bloch sphere) of a two-level system 
driven by Eq.JS)) in a time interval t G [0,4]. Both R and v 
range from to 7. The axes Px,Py,Pz are the Bloch vectors. 
The initial state is |e), i.e., pz = l,Pcc = Py — 0, and — 
= 0. Note that there is a set of unaccessible states on the 
other side of the Bloch sphere, which can not be seen from 
this angle, (a) C = 2, (b) C = 6, (c)C = 8, (d) C = 10, 
(e)C - 12, (f)C = 14, (g)C = 20 and (h)C = 100. 



By the unitary transformation F(^, 0), the Hamilto- 
nian Hni is transformed into. 
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(4) 
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+ (^ cos^ + Rni COS sin cos 4* ~ 2 ^^^^ ^ '^4)^x 



{Rni COS sin ^ sin — - sin^ sin 2(j))cry , 



(5) 
(6) 



where Rni = R — C{ilj\az\i^). We have performed exten- 
sive numerical simulations for the Schrodinger equation 
ih-^l^p) = H'^j\ip) with ^ = = 0, select results are pre- 
sented in Fig|2l Two observations can be made from Fig. 
[21 (1) The nonlinearity affects the controllability of the 
two-level system, regardless of how small the nonlinear 
coupling constant C is, (2) the larger the nonlinear term 
is, the smaller the set of the accessible states. Further 
numerical simulation shows that 9 and (j) (determining 
the initial state) can change the accessible set of states 
but not the two observations, this is shown in Fig|3l 

The nonlinearity in nonlinear quantum system may 
come from feedback control ^ {^1)\cf z\'4^) cr z - This feedback 
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FIG. 3: (color online) The same as Fig|2] but with — 
7r/4,(/) = 0. C — 20 was chosen for this plot, (a) is the same 
as (b), except the position of axes px and py. 



(c) 




control may be understood as applying a pulse to shift 
the energy level of the quantum system, the amplitude of 
the pulse is proportional to a detection event with a rate 
C described by {iIj\cf z\il^) - To overcome the uncontrolla- 
bility induced by the feedback, we introduce a Lyapunov 
control f{t) to replace the coupling constant C in the 
nonlinear term z\i^) (J z - The Lyapunov control may 

be designed according to the Lyapunov control theory as 
follows [1^. For a two- level system, its Hamiltonian can 
be rewritten as. 



Hni = Ho f{t)Hi, 



(7) 



where Hq is the free evolution Hamiltonian and Hi is the 
control Hamiltonian. The general control task we con- 
sider can be formulated as, given a target state \tljd{t))^ 
we wish to apply a certain control field f{t) to the system 
that modifies its dynamics such that l'^d(^)) as 

t ^ oo. Since the free Hamiltonian can in general not 
be turned off, it is natural to assume \ipd{t)) to be time- 
dependent and satisfies 



in^^\Mt)) = ifo\Mt))- 



(8) 



Since the evolution of both \ip{t)) and \ipd{t)) are unitary 
in our case, we can define a function 



v[\Mt)),\m)] = i-\{Mt)\m)\' 



(9) 



to measure the distance between the resulting and target 
states. Clearly V > with equality only if \il^d{t)) = 
\ip{t)). Taking derivative of V with respect to time t, we 
have {\^l^d) = li^dit)) and = \ip{t)) hereafter) 



V 



-2f{t)lm{{^d\Hr\^){^m). 



(10) 



where Im(...) denote the imaginary part of (...). So when 
we choose f{t) = Klm{{tljd\Hi\tlj){ip\ipd)) with a rate hz > 
0, we have V < 0. Therefore F is a Lyapunov function 
for the following dynamical system, 

d 

fit) = Klm{{i,d\Hi\i;){i^\i^d)). (11) 



FIG. 4: (color online) Accessible states represented by the 
Bloch vector p = {Px,Py,Pz)- The two- level system is driven 
by Ea. (|13p with Lyapunov feedback control /(t) evolving in 
a time interval t G [0,4] with — (j) — 0. R and v range from 
to 0.5. (a) At = 0, (b) = 3, (c) At = 9, (d) At = 27. 



For a two- level system, \ipd) always can be written as 
\ijd) = c{t)\e) + d{t)\g) with |c(t)|2 + \d{t)\^ = 1, and 
IV^) = a{t)\e) + b{t)\g) with \a{t)\^ + \b{t)\^ = L When 
the control Hamiltonian takes Hi = it is easy to 

find that the Lyapunov control f{t) = — 2/^Im(c(i*a*6). In 
the following, we shall focus on the control Hamiltonian 
Hi = ^(V^lcr^lV^)^^, which yields the Lyapunov control 



f{t) = -/^mlm[(ac* + 6d*)(a*c-6*d)], 
m = |ap-|6p. 



(12) 



where we omitted the argument t of a{t) ^ b{t) ^ c{t) and 
d{t) to shorten the notations. Clearly, the Lyapunov con- 
trol renders the dynamics of the quantum system nonlin- 
ear even if the control Hamiltonian is linear. We have 
performed extensive numerical simulations for the dy- 
namics of these nonlinear system, the numerical simula- 
tions show that the two- level system described by Eq.([T]) 
with Lyapunov control f{t) is controllable, namely an 
arbitrary pure state is accessible driven by 



Hni 



R 



(13) 



Then a natural question arises for how does the rate 
K in Lyapunov control f(t) affect the accessible set of 
states? Fig. |4] shows the sets of accessible state reached 
by controlling the parameter R and in a small range. 
In other words, the control parameters R and v are re- 
stricted in a regime much smaller than that by which the 
two- level system is controllable. We find that the rate 
affects the accessible states. For n below a critical value 
/^c, the larger the n is, the bigger the set of accessible 
states. For n > the situation changes, smaller k fa- 
vors the set of accessible states. In FiglSl we plot the 
short-time behavior of the accessible states, observations 
similar to FigHJcan be found. A common feature we find 
from Fig m and [5] is that as n increases, states near the ini- 
tial state |e) become easy to access. This finding depends 
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FIG. 5: (color online) The same as FiglH but the time interval 
is [0, 0.3], and R and v vary in [0, 7]. (a) At = 0, (b) = 9, (c) 
K = 81, (d) K = 243. 

on and (f). Finally, we address the convergence for the 
control system. By the Lasalle's invariance principle [2 Ij, 
the largest invariant set is empty, so there is not any in- 
variant set for the problem under consideration. The rea- 



son is that we choose \iljd{t = 0)) = = 0)) throughout 
this paper. 

In summary, we have investigated the nonlinear effect 
on the controllability of a two-level system. This non- 
linear effect can turn a controllable quantum system un- 
controllable. The accessible set of states under nonlinear 
effect depends on the nonlinear coefficient C and the ini- 
tial state. To overcome this uncontrollability induced by 
the nonlinear effect, we propose to use Lyapunov control 
to manipulate the two-level system, Lyapunov function 
for the control system is constructed and the dependence 
of accessible set of states, which can be reached in a short- 
time limit and within a small range of control parame- 
ters on the rate are shown and discussed. This study 
suggests that feedback control that can induce nonlin- 
ear effect changes the controllability of quantum system, 
Lyapunov control is better in this case for manipulating 
a quantum system. 

This work was supported by NSF of China under 
Grant No. 10775023. 
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